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NEW ASPECTS OF SUBSONIC AERODYNAMIC NOISE THEORY 
by Marvin E. Goldstein and Walton L. Howes 
Lewis Research Center 

SUMMARY 

A theory of aerodynamic noise is presented which differs from Lighthill's theory pri- 
marily in the way in which convection of the noise sources is treated. The sound direc- 
tivity pattern obtained from the present theory agrees better with jet -noise directivity 
data than does that obtained from Lighthill’s theory. The resxalts imply that the shear - 
noise contribution to jet noise is smaller than previously expected. 


INTRODUCTION 

The theory of aerodynamic noise was developed by Lighthill in two classic papers 
(refs. 1 and 2). Lighthill showed that, in the absence of solid boundaries, aerodynamic 
noise is produced in a fluid by turbvilent stress fluctuations acting as acoustic quadru- 
poles. In a moving fluid these quadrupoles are carried along by the mean flow at some 
convection Mach number M . For nearly parallel flows (as in jets), convection of the 
quadrupole sources introduces a factor (1 - M cos 0 )~ into the noise directivity pattern, 
where 0 is the angle measured from the mean flow direction. In Lighthill’s theory this 
convective directivity factor dominates the directivity pattern of jet noise. However, it 
has been found that, when experimentally measured values of M are used in the theory, 
the predicted directivity patterns do not give good agreement with experiment. 

In the present report a possible alternative theory of jet noise is proposed. This 
theory leads to a formula for predicting noise which differs from the formulas obtained 
from Lighthill’s theory in that the exponent -5 in the convective factor is replaced by the 
exponent -3. This gives better agreement with existing experimental directivity data. 

Since the present theory begins with a convected wave equation, it starts from a 
more general point of view than Lighthill’s theory. The effects of convection and refrac- 
tion of the emitte d sound waves by the mean flow are neglected, and the resulting wave 

^he exponent -5 was first obtained by Ffowcs Williams (ref. 3), Lighthill errone- 
ously obtained an exponent -6. 



equation is then solved for the far -field sound pressure by using the free -space Green’s 
function. Nejct, the source term (which involves the fluctuating Reynolds stress in es- 
sentially the same way as Lighthill’s source term) is evaluated by using the isotropic 
turbulence model adopted by Lilley (ref. 4) and by Ribner (ref. 5). Finally, it is shown 
that the directivity patterns of the far -field sound intensity obtained from this theory are 
in better agreement with experimental data than those obtained from Lighthill’s theory. 


DERIVATION OF BASIC EQUATION 


We begin by deriving an approximate wave equation very similar to Phillips exact 
wave equation (ref. 6). However, it is convenient to derive this equation from first prin- 
ciples, instead of starting with Phillips equation, since certain approximations which we 
shall make are most conveniently introduced durir^ the derivation. As in reference 6, 
we start from the Navier -Stokes equation 


DV 

Dt 


i V p + i V . e 
P P 


and the continuity equation 


V • V = 0 

p Dt 


where 


. ( 1 ) 


\ 


( 2 ) 


Dt 


-L+ V • V 
0t 


p is the pressvire, p is the density, V is the velocity, and e is the viscous stress ten 
sor. (All symbols are defined in appendix B. ) We also assume that the ideal -gas 
relation 


p = Rp © (3) 

holds. The velocity V can be expressed as the sum of a time -averaged part U and a 
fluctuating part v’. Thus, 


V = U + V’ 


(4) 
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Furthermore, suppose that the mean velocity is parallel to the Xj -coordinate axis and 
varies only in the Xg-direction in a locally Cartesian coordinate system (xj, Xg, Xg). 
Since refraction effects will be neglected, this is a reasonable assumption for a jet. 
Thus, 


U = iUCxg) 


( 5 ) 


The density can be expressed as 


p = Pq + p’ (6) 

where p^ is the time-averaged density and p’ is the fluctuating part of the density. Al- 
though we shall allow the maximum mean-flow Mach number to be any value less than 
unity, the Mach number of the turbulence will always be much less than unity. In fact, 
for jet flows, the Mach number of the turbulence will be less than 0. 1 or 0, 2 (ref. 7). 
Hence, the tvirbulent flow is essentially incompressible. And the fluctuating part of the 
density is essentially associated with acoustic waves and not with the turbulence. 

As pointed out in reference 8, the interaction of the emitted sound with the turbu- 
lence is negligible. This interaction is accounted for in the equations by the terms which 
involve products of the turbulent quantities with acoustic quantities. Such terms will, 
therefore, be neglected. We shall also neglect any terms which involve squares of acous- 
tics quantities since these terms are assumed to be small. 

It follows from equations (4) to (6) that the exact time-averaged form of the continu- 
ity equation (2) is 


Since V' is the sum of turbulent and acoustic velocities, the assumptions discussed in 
the preceding paragraphs imply that the right side of this equation can be neglected. 
Therefore, assuming an initial uniform density upstream, the time -averaged density p^ 
will be independent of position. 

It also follows from equations (4) to (6) that 

51£ = + U ^ + V’ • Vp’ 

Dt at 3Xj 


The assumptions discussed in the preceding paragraphs again imply that the last term on 
the right side of this equation can be neglected. We can therefore write 
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Dt at dx. 


so that the continuity equation becomes 

1 Ie + u -^j+ V . V = 0 
p at axw 

Now it can easily be shown from the formula for entropy change in an ideal gas 
(ref. 9) that 


^ = JLdp - -LdS 

p yp Cp 

where S is the entropy, Cp is the specific heat at constant pressure, and y is the 
ratio of specific heats. Hence, 

V • V =-/"— + U — Y 
yat ax J 



(7) 


where we have pvit 


n s 1 In P- (8) 

r Po 

and have, for convenience, chosen p to be the average far -field pressure. The flue- 
-► o 2 

tuating velocity V can be decomposed into a turbvilent part w plus an acoustic part 

V, where v is called the acoustic particle velocity. Thus, 

V' = w + V (9) 

Since the turbulence is essentially incompressible, we can suppose that 

V • w = 0 (10) 


Hence, 


2 


Notice that we are ignoring acoustic streaming. 
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V-V = V.V' = V.V 


Since v is an acoustic quantity, this shows that V • V is essentially an acoustic quan- 
tity. In view of the assumptions previously described, we can, therefore, write 


D 

Dt 


V • 



V 


Hence, it follows from equation (7) that 


^ V • V = 
Dt 


^ + u -L-f n + J_/i. + u -LY s 


I0t 3x 


1 


Sr 


Upon using the vector identity 


V • V = 
Dt 


9 DV. 0V. 0V. 


0Xj Dt 0x. -0Xj 


given in reference 6, this becomes 




0V. 0V. 
0x, 




where the summation convention has been adopted. After substituting equations (4), (5), 
and (9) into this equation and neglecting the squares of acoustic quantities and the prod- 
ucts of acoustic and timbulent quantities, we obtain 


0 


DV. 


0x. Dt 




S + 2 


dU 

dx„ 



0W. 0W. 
0Xj 0Xj 


( 11 ) 


where we have put 


Wi = 


6iiU+Wi 


( 12 ) 


with 6jj the Kronecker delta. Furthermore, using the fact that, for an ideal gas, the 
speed of sound c is given by 


5 



the Navier -Stokes equation becomes 


c2=yP 

P 


DV 

Dt 


= - c^vn + i V 

P 


( 13 ) 


Hence, upon taking the divergence of equation (13) and using the result to eliminate the 
left side of equation (11), we obtain 




n - V • c^vn 


dU ^2 
dx^ ^ 


3Wj 0Wj 

0x. 3Xj 




+ -i- /L + U S 
Cp \^3t dx^j 


(14) 


This equation differs from Phillips' equation, (eq. (2. 8) in ref. 6), which is exact, in two 
important respects, namely the operator D/Dt in Phillips' equation is replaced by the 
operator 


2_ 

3t 


+ 



and the terms 






in Phillips' equation are replaced by the terms 



dxg 3Xj 


3W. 3W. 
3x^ 3Xj 


The derivation carried out above shows that these approximations amount to neglecting 
squares of acoustic qviantities and the interaction between acoustic waves and turbulence. 
As pointed out in reference 8, for subsonic flow, the neglect of the interaction of the 
sound with the turbulence is justified on the basis that the acoustic waves emitted by the 
turbulence are of such long wavelengths compared with the size of the turbulent eddies 
that they pass right through the turbulence without being significantly affected. 


6 



The changes in entropy are related to the heat conduction through the energy equa- 
tion. However, for the high Reynolds numbers which are of interest in jet flows, the di- 
rect effects of heat conduction and fluid viscosity are likely to be unimportant for cold 
jets. Hence, we shall assume that the entropy is constant and that the speed of sound has 
the constant value c^. In addition, we shall neglect the viscous stress tensor e. Then, 
equation (14) becomes 


fl + u -i- 

,9t 3Xi 


n - CqV^h - 2 


dx2 0Xj 3Xj 0Xj 


(15) 


Upon substituting equations (4), (9), and (12) into equation (13), neglecting terms involv- 
ii^ squares of acoustic quantities and interaction of the sound with the turbulence, the 
Xg -component of equation (13) becomes 


(l.U-l 

l0t 


V2 = -c2iIL- l.W.^)W2 

0Xj 1 ^ ° axg \0t ^ 0xj ' 


(16) 


Then, by eliminating Vg between equations (15) and (16), we obtain 


(1 + U-l. 
l0t 3Xi 


(1 + U- 
\0t 0X, 


n - CqV^h 


+ 2c 


2 dU 0^n 

dXg 0Xj0X2 




0W, 


+ W. 


)0t 


dx2 0Xj \ 0t 


0Xi 


(17) 


Now, it follows, from equation (12) and the fact that U is a function of X 2 only, that the 
right side of equation ( 17) can be written as 


(0t 0XjJl 0Xj 0x. 


+ 2 


dU ^^2 ^ 
dx^ 0x 


2 ^ dll I 2 ^ TJ 
0Xi 0*2 I 0t 0Xj 


+ W. 


0W, 


0X. 


^ + U 




^0t 0Xjyy0Xj 0Xj 


dx2 0Xj 


0Xi^ 
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However, it follows from equation (10) that 


w. 


Bwg 9 (wjW2) 
9x^ 3x. 


and 


3Wj 0Wj 0 w.Wj 

0X. 0X. 0X.0X. 

J * J 


Hence, equation (17) becomes 

V2 


1 + U-l 

l0t 0X 


ly 


^ + u -2-1 n - c“v^n 

l0t 0X, 


+ 2c 


2 dU 0^n 
dx2 0x^0X2 


= ^L+U 


y0t ^1/L j 


0^(WiWj) 


dU 


dXr 


Since 




0X 


J 0Xj0Xj 0Xj |0Xj J ^2 


U 


3(WjW.) 


dU SCwgWi) 


0 

0X, 


^ (Uw.w.) -2 ^^(w.W 2) -w| 


2 d^U 


0Xj^03^ 


dx 2 0x. 


dx' 


2J 


and 


i_ uv^n = -L 


0X, 


0X, 


vun-2^iIL-^n 

*^2 ^^2 dx| 


this equation can be written as 


0^(w2Wj) 

0Xj0Xj 
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where we have put 



ll\[rERPRETATION OF BASIC EQUATION 

Equation (18) is the starting point of the present analysis. We shall suppose that the 
turbulent and mean flows are given. Then, since the terms on the right vanish outside of 
the jet, we can treat the right side of equation (18) as a soxirce term. Since no inter- 
action between the sound wave and the turbulence is allowed once the sound is emitted, 
from the acoustical viewpoint each turbulent eddy acts independently. The sound is sim- 
ply convected and refracted by the mean flow. In fact, the left side of equation (18) is a 
wave equation for the quantity r in a parallel flow with velocity U varying only in the 
Xg -direction and with the term accovmting for the direct refraction of the sovind by the ve- 
locity gradients neglected (see ref. 10). The source terms on the right side represent 
the production of sound by the turbvilence. The consequences of equation (18) coincide 
with those of Lighthill’s theory if the mean flow is negligible, but otherwise (as we shall 
see) equation (18) can lead to different results. 

NEGLECT OF CONVECTION AND REFRACTION EFFECTS 

The operator ^9/3t + U d/d'x^^ on the left side of equation (18) is the time derivative 
in a coordinate system moving with the local mean velocity U, which we shall see is 
roughly equal to the velocity at which the individual sound sources move. This time de- 
rivative therefore roughly has the effect of multiplying the term on which it operates by 
the average angular frequency Q of the soxmd in a coordinate system moving with the 
source. Hence, if the wavelength 2irc^/^ is large compared with the jet diameter, the 
first term on the left side of eqxiation (18) can be neglected by comparison with the second 
term at all points within the jet. On the other hand, since U is only nonzero within the 
jet, the operator in the first term reduces to 3 /9t outside the jet. However, the 


9 



operator 3/3t roughly has the effect of multiplying the term on which it operates by the 
average frequency u> of the sound in a fixed frame. Hence, if the typical wavelengths 
2irc^/(ii and 2irc^/Q, are large compared to the jet diameter, we can replace the oper- 
ator on the left side of equation (18) by the free -space wave operator. 


Now the data in references 11 and 12 show that these typical wavelengths vary from 
around 6 to 10 jet diameters. Hence, we shall make the assumption described in the pre- 
ceding paragraph and replace equation (18) by 


^l-cVir = 

^3t^ 


^0x^9x. \9t 


- + u j_| (W.W.) 


- 4 


0X, 




d^U 


LdXg 


( 20 ) 


This approximation amounts to neglecting the convection and refraction of the sound by 
the mean flow; an approximation which is also made in Lighthill's theory. 


SIMPLIFICATION OF SOURCE TERM 

Within the jet the acoustic velocity is certainly negligible compared with the flow 
velocity 


W. = 6iiU + Wi 


( 12 ) 


Thus, within the jet the X 2 -component of the momentum eqmtion (16) can be written as 

And since the source term in equation (20) vanishes outside the jet, we can substitute this 
expression into the source term to obtain 
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( 21 ) 


^--cy\v = 

V3t2 


+ U -i-\(w.w.) 


dx^dx^ \ dt 3Xj 


1.U-LK--L 

9Xj dX2 y9t ^^1/ 


d^U / 2 2n\ 

72 i"2 ^ 




We now introduce the vector potential A by 


I(J,t)=lvx I^Miy 
ij J \x. y| 


( 22 ) 


Then, as shown in books on electrodynamics (e.g. , ref. 13, sec. 5.3), since V • w = 0, 
this equation implies that 


w = V X A 


(23) 


and 


V • A = 0 


(24) 


Thus, in particular , 


g 3Ai 



where is the permutation symbol. Upon substituting this into eqviation (21) and re- 

calling that U is a fxmction of Xg only, we obtain 


/ii - c vV = 

A + U ^ \/w w + 46 e , A, Y 

_ _3_ 

^ fw5 + c^n^ 


[ax,axj^at >=jj 

0Xj 

j 


Now, as suggested by Lilley (ref. 4) and demonstrated by others (refs. 14 and 15), 
in a jet shear layer the mean velocity gradient dU/dXg (mean shear) is slowly varying 
over the narrow strip along the center of the mixing region where most of the turbulent 
energy occurs. Hence, terms corresponding to the variation in dU/dXg should only 
make a small contribution to the source term and, therefore, will be neglected. Thus, 
we will assume that 


11 



= constant 


dU 

dXg 


Then equation (25) becomes 


° ) 





where we have pxit 




( 26 ) 


(27) 


This equation is similar to Lighthill’s equation in that it is an ordinary inhomogeneous 
wave equation for a viniform medium containing a quadrupole soiirce. However, it differs 
jfrom Lighthill’s equation in the important respects that the dependent variable is T in- 
stead of the condensation (i. e. , (p - p^)/p^) and the quadrupole source term is consider- 
ably different from Lighthill’s. 


RETARDED POTENTIAL SOLUTION OF BASIC EQUATION 


The retarded potential solution of equation (26) for r which vanishes at infinity is 


r(x) = 



-i+U(y2)-^ 

9t ^ 6yj 



(28) 


where r= |x - y( and the vectors y and x now denote sovirce and receiver points, re- 
spectively, In addition, S/Sy^^ indicates partial differentiation with respect to y^ with r 
held fixed. Now, for any function F(y,r,t), 


6F _ 3F 9F 3r _ 9F ^ 8F 

6yj 9yj 9r 3y. 3y. 3Xj 


where d/dy^ denotes partial differentiation with x fixed and 0/3x, denotes partial dif- 
ferentiation with y fixed. Then, 
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6^f e^F ^ a^F ^ ^ a^F 

SYiSy. 3y.0yj 0y.0x. 0y.0x. 0x.0x^ 


Upon inserting these results in equation (22), using the divergence theorem and the fact 
that 


we obtain 


Since 


^ik 




as 



CO 





dy 



- y? 



it follows from eqviation ( 19) (since U vanishes outside the jet) that, when r is large 
enough so that x is in the radiation field of each turbulent eddy (so that terms of higher 
order in 1/r can be neglected, see ref. 1), the preceding equation for r reduces to 


* J 


U(y 2 ) • (x - y) 

V 




t 



(29) 


Since (p - Pq)/Pq is certainly very small in the radiation field, it follows from equa- 
tion (8) that 


n 


P-Pq 

Po^'o 


which is eqxial to the condensation in the radiation field. 


(30) 
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SOLUTION IN TERMS OF SPECTRA 


We shall assume that the turbulence is statistically stationary in time. Then we 
shall suppose (as usual, in order to ensure the convergence of the Fourier integrals 
which will be introduced below) that = 0 for jt| > T, where T is some large time 
which will be put equal to infinity at the end of the analysis. We now introduce the 
Fovirier transforms 






e Ty dt 


Then, upon using the asymptotic approximation (30), the Fourier transform of equa- 
tion (29) becomes 

4i7Cg J r® 


1 - 


U(y 2 ) • (x - y) 




Tij dy (31) 


In the far field of the entire flow, the spectral density of the intensity is given by 
(ref. 16) 


I 

O) 


(x) = 



Substituting equation (31) into this relation and using the convolution theorem gives, in 
view of equation (27), 



^ U<y2)(x-y) 

U(z2)*(x-z) 

32^5 JJJ 

<=ol*-y| 



7)dy dz dr 


(31a) 
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where the operation of taking the real part is certainly justified since is real. (This 
operation is performed for convenience in subsequent manipulations.) Also, 

^ = w.w.w^w^ - W.W. w^w’ 


+ 166 jj 6 jj^€ 2 jm^ 2 Zn 





dU 

dYo 




(32) 


is a fourth -order, two -point, time -delayed barbulence correlation tensor where 
7 ) =z ~y. The primes indicate that the quantities are to be evaluated at z = y + ^ and 
t + T, while the unprimed quantities are evaluated at y and t. The overbar denotes the 
average 


f(t,t + t) = 

2T 



f(t,t + 7)dt 


The term w.w. w[.w’, , which does not contribute to I, ,(x), has been introduced for con- 
venience. Upon changing the variables of integration from y and z to y and ^ and 
collecting terms, the equation for I^(x) becomes 




4 

Po 





U(y«) • (x-y) 
1 - ^ - 

U(y2+T)2) ' (x -y -’Ji 

JJJ t?-yl^|x-y -nP 

Co|x-yj _ 

cJx-7-n] 




(33) 


FAR-FIELD EXPANSION OF SOLUTION 


Because the distance | rjj over which the correlation is nonzero is certainly 

smaller than any overall dimension of the region of turbulence, we can always suppose 
that the observation point x is svifficiently far away from the flow so that 

|x - y| » (^1 


Therefore, vq)on expanding the integrand in equation (33) and neglecting terms of order 
j^j/r, we obtain 


fff W 

JJJ r® 


e 

6 


, "(yg) 

1 cos e 


^ U(y2_^) 


cos 6 




(34) 
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where we have put 


r s X - y 

ri = X, - Vi 


and 


r rU 

is the cosine of the angle 9 between the direction of mean flow and the direction of ob- 
servation. 

Now, 


U(72 + V2^ - ^(y2^ ~ — ^2 ”^ ! '. i •• 

dY2 

since dU/dyg is assumed to be constant. For the range of |tfj over which the correla- 
tion is nonzero, the ratio of the second term to the first term on the right side of 

this expression is of the order of the ratio of the transverse correlation length to the 
thickness of the mixii^ region, since for jets the change in 11 ( 72 ) occurs over the width 
of the mixing region. It can be seen from the resvilts presented in references 14, 15, 
and 17 (also see discussion in ref. 7) that this ratio is of the order of 1/6 to 1/10. 

Hence, we shall neglect the second term compared with the first and replace U(y 2 + 
by U(Y 2 ^ equation (27). Since this U now represents a sort of average velocity over 
the eddy, we shall set it equal to the "eddy convection velocity" Ug(y 2 )- fact, as 
pointed out by Lilley (ref. 4) and verified experimentally in reference 15, the chief noise- 
emitting eddies are confined to a narrow strip in the center of the jet mixir^ region; and 
in this strip (ref. 15) the eddy convection velocity is very nearly equal to the local mean 
speed U. (See refs. 15 and 17 for a discussion of the convection velocity.) Upon collect- 
ing terms, equation (34) now becomes 


frr r^r.ri^r; iw[(r.7]/rc )-t] 

[^(J)= ° /?«/// ^ ^ ^ ^ e ° (l-M^ cos 0 r (y, 77 , T)dy d^ d t 

327r^c® JJJ r® 

° (35) 


16 



where 




is the eddy convection Mach mrniber. 


INTRODUCTION OF MOVING COORDINATES AND NEGLECT OF RETARDED TIME 


When Tj is transverse to the flow, the changes in r • rj/rc^ over the "eddy" cor- 
relation distance are small compared with the change in t over the eddy decay time. 
Hence, we can neglect the term r • if/rc^ in 


ejq) 




However, when tj is in the direction of the mean flow (see ref. 16), the changes of the 
correlation with respect to time are much more rapid since the eddies are moving with 
the flow. In order to compensate for this, we introduce a set of coordinates which move 
with the convection velocity U^. Thus, we introduce the new variable 

?=^F-UJy2)T (36) 

where is in the yj-direction. Following Ff owes Williams, we define the moving - 
axis correlation function by 


(37) 

Upon introducing this change of variable into the integral in equation (35) and noting that 
the Jacobian of the transformation is unity, we find that 


4 




cos e] 


2 f -iw(l-M cos e)T f (ioj/c r)T-r 

fl^ *= /e ° Rrjki(y,«,7)d«dTdy 


(38) 
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Now let I denote a typical correlation length of the turbulence, let Mg be the char- 
acteristic Mach n;imber of the tvirbulence, and let denote a typical decay time of the 
moving -axis correlation. It is shown in reference 7 (see sec. 3, eq. (17)) that 


— «M^ 

■^f^o 

However, for jet flows the maximum turbulent velocity is about one -tenth, and certainly 
not more than two-tenths, of the maximum mean-flow velocity. Consequently, for sub- 
sonic flows 


This shows that | • r/c^r is negligible in comparison with (1 - M^ cos 9)r. Hence, we 
can set | • r/c^r equal to zero to obtain 


= 


4 

^ PQ 


fie 


f 

J 




Mj, cos e 


vp f* -iw(l-M cos B)t r ^ 


dTdy 


(39) 


Following Ribner (ref. 5) we let I^(x|y) denote the intensity at the point x caused by 
the sound emitted from a unit volume at y. Then, 

and it follows from equation (39) that 




r.r.rj^r^ 
32jt^c^ r® 





( 40 ) 


where 


= o)(l - Mg cos 0 ) 
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is the frequency which would be observed in a frame of reference moving with the eddy 
convection Mach number M^. 

We can now write, in view of equations (32) and (37), 


'dU' 


' 2y 


dU 


(41) 

2 


where 


r) - Wjw.w^w’ - WjW. w^w’ 


Q„,,kz(y,?,^) ^AmW]^w’ 


0^ (y,^,r)-A^A; 


(42) 

(43) 

(44) 


We shall now show, under the relatively mild restriction that the turbulence is 
locally homogeneous and incompressible, that the volume integrals of the correlations 
and can be expressed in a simple way entirely in terms of two -point ve- 

locity correlations. 


LOCALLY HOMOGENEOUS TURBULENCE; REDUCTION 
OF CORRELATION VOLUMES 

The theory of incompressible, locally homogeneous turbulence is developed in ref- 
erence 18 and in chapter 4 of reference 19. The assumption of local homogeneity implies 
that any two-point, two-time correlation function, say Q(y,z, t), is a function only of 
7j = z - y and t. 
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It follows from equation (22) that 








where the doiible primes denote quantities evaluated at y" and t and where 
7 ^^^ = z - y”. Or, since the Jacobian of the transformation [where 

- Uj.(y 2 )T] is unity, 


^,kZ ■ 


L 'mrp y ®p, 




where 




is the two-point, time-delayed, triple -velocity correlation. Thus, 


^ %i,U " 4 ^ ^mrp 




94^^ 


i di" d4^^ 


Hence, upon using equation (Al) derived in appendix A, this becomes 


f 0« = J'n.rpy 4” \,U <4“’ ' 


«rRp, (<« 


It is shown in reference 20 that this integral is indeed convergent (see remarks immedi- 
ately proceeding equation (6. 17) of ref. 20). 

In order to evaluate the volume integral of notice that it follows from equa- 

tion (22) that 


20 



1 — JT _ ^mpq^nrs 
167T^ 


^mpq^nrs 






|y - y"l~^ — ^ |z - dy" dy”* 


ayi.” 


|y - y"(-l _i_ |z - y” - dy” 


dTj. 


( 2 ) 




where the triple -primed quantities are to be evaluated at y'" and t + t and where 
jj(2) = yttt _ y»'. The integral in square brackets was evaluated in the appendix of refer- 
ence 21 and found to be 


27T 





Hence, 


a„,a; 

m n 


mpq^nrs 

87T 


w”w"’ 
q s 




V' 





Or since the Jacobian of the transformation —77^^^ j^where 5 77^^^ - Ug(y2)'^J 
imity, 


r\ - ^mpq^nrs 
8v 


/ 


R 

qs^^ > 


3412)34(2) 


d|(2) 


where 


Rij(4 , t) = w^wJ 


( 47 ) 


is the two -point, time-delayed, second-order velocity correlation tensor. Hence, 
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The doxible integral was evaluated in appendix I of reference 4 by using the methods 
of reference 21 and found to be 


— 

30 








Hence, 



(F, T)d? = 


^mpq^nrs 

30 




2^p?r ^pr 






(48) 


It is shown in reference 20 that this integral is indeed convergent (see remarks on top of 
p. 386, ref. 20). 

Upon using the identity and the fact that = -e^kj, 

it follows from equations (41), (46), and (48) that 




8 


m ^ ^li^lk(%Sjq - 62qy(62r5^s 


■ ^2shr^ 


^dU> 

&2j 


Ji^Vv * - I hi ^ f(h\u - h^2, 


(49) 


The first two integrals on the right side of equation (49) represent the self noise and the 
shear noise, respectively (refs. 22 and 4). This terminology was introduced by Lilley in 
reference 4 to indicate that the former term represents noise generated by turbulent - 
turbulent interactions, whereas the latter term represents noise generated by turbulent - 
mean-shear interactions. The last integral represents a coupling between the shear 
noise and the self noise. 

Equation (49) cannot be further simplified without introducing some additional as- 
sumptions about the turbulence. We first introduce the joint normality approximation. 
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REPRESENTATION OF FOURTH-ORDER CORRELATIONS 
IN TERMS OF SECOND-ORDER CORRELATIONS 


It is argued by Batchelor (ref, 18) that the part of the joint probability distribution of 
the velocities at a fixed time associated with the energy -bearing eddies is approximately 
normal, at least insofar as the velocities at two points are concerned. This approxima- 
tion is better for some purposes than others. Thus, in reference 18 (p. 176) it gives 
reasonably accimate predictions about the relation between the second- and fovmth-order 
correlations. This relation (ref. 18, eq. (8. 3. 11) is found to be (see ref. 23, sec. 2.1.7 
for derivation) 


WiWjW^W^ = ^i^j 


W^W^ + W.WiV WjW', + w.w', w^w,'. 




at T = 0 


(50) 


But, by extending the reasoning used by Batchelor in section 8.2, we can argue that, 
when the velocity correlations in the moving frame are separated in time as well as in 
space, their correlation will be subject to even more random influences from the neigh- 
boring flow than when they occur at the same time. In accordance with the central limit 
theorem , these influences will tend to ensure fiirther the normality of the joint probability 
distribution. Hence, we have even more reason to expect equation (50) to be valid with 
T * 0 , and we can now write from equation (42) 

^ijkz = ^ik^y ’ 

so that equation (49) now becomes 




8 


ik^jZ ^ ^iZ^jk^^^^ - fg Sli«lk(S2pSjq - ^2q«jp^(52rSzs 


^2s®Zr^ 


^ dU 


f{%ir * V«>qs '>« " I ^11 5 


(51) 


RESULTS FOR ISOTROPIC TURBULENCE MODEL 

In order to deduce the jet -noise directivity pattern from equations (40) and (51), it is 
necessary to make some additional assumption about the turbulence. Perhaps the sim- 
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plest plausible assumption is that made by Ribner and Lilley (refs. 4, 5, and 22) namely 
that the turbulence is locally homogeneous and isotropic. A more refined model for the 
turbulence (assumption of axisymmetric turbulence) is given in reference 24 . However , 
since (as we shall see) the convective factor dominates the overall directivity pattern 
(i.e. , at least insofar as the overall intensity itself, but not its spectral density, is con- 
cerned), it is felt that the assumption of isotropic turbulence will suffice for predicting 
the directivity of the intensity. We shall, therefore, assume that the moving -axis corre- 
lation tensors are isotropic tensors. Thus, it is shown in reference 18 (p. 42) that the 
general third -order isotropic tensor which is symmetric in its second two indices is 

where C, D, and E are fvmctions of y and t and are fimctions only of the magnitude 
^ of J and not its direction. Hence,, we put 


Then 


^2^j,kZ ■ ^i^2,kZ ■ ^^^2^k®Zj ' ^j^k^Z2 " ^Z^2^ 


And, therefore, upon noting that only even functions of can contribute to the integral, 
we find that 






+ ^Z2 V 




(no sum on j) 


= 0 (52) 

since ^ d^. This shows that the coupling term in equation (51) 
vanishes. 

It is shown in reference 24 that the second -order isotropic correlation tensor is 

Hij(y,?,T)=A4.Cj+B6.j (53) 

where A and B are functions of y and t and are functions only of the magnitude | 
of 4 and not of its direction. It is also shown in reference 24 that the requirements of 
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continuity dictate that the functions A and B are related to a single function F by 


A - 1 1 9F ^ 

Mi? 

B = Fa8I 
2 0?^ 


( 54 ) 


Hence, upon noting that only even functions of contribute to the integral, we find that 
(«2p«iq - «2q«ipX«2r«is ‘ «2s«ir' / (^Spir ^ 


= 26q 


n 


f H * «>22 -f * h/hi 


(no sum on j , Z ) 


26j;(l 


- V/[(^^^ - 2«ll2Rl2 


d^ (no sum on j) 


= ?. 6j^(l - Sgj) 


where we have used the fact that for any function f(^), 

J f(?)|2 J d? 

Hence, upon integrating by parts, we find that 

^^2p®jq “ ^2q®jp^^^2r^Zs " J (^^p^r ■*■ ^pr^ )^qs 


= - J - «2J> / - «2i> / «2"*11 d^ (no sum on j) (55) 
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Finally, it follows from equations (53) and (54), after applying some routine manipula- 
tions given in appendix C of reference 25, that 

Hence, upon substituting equations (52), (55), and (56) into equation (51) and usii^ the re- 
sult in equation (40), we obtain 


r.2 2 

— ^ n Cl) 

— 5T7 




R?, dl'dr - 1® cos^ e (cos^ e + sin^ 0 sln^ 
“ 3 




dT 


(57) 


where the azimuthal angle cp is given by sin e sin (p = r^/r. 

EQUATIONS FOR AXI SYMMETRIC FLOWS AND TOTAL INTENSITY 

The total intensity I(x|y) at each point is obtained by integrating the spectral density 
of the intensity over all frequencies. Thus, 

I(x|y) = /; 


Substituting equation (57) into this expression shows that 


I(x|y) = 


1677^c^r^ (1 - Mp cos 9)^ 


' f 

_ 


Rjl dM - — cos^ 0 (cos^ 0 + sin^ 0 sin^ (p) i- 

4=0 \ 


duV 

J 


^ 2^11 


^7=0J 
(58) 


As indicated by Ribner (refs. 5 and 22), for round jets it is appropriate to average the 
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intensity over the azimuthal angle (p. When this is done, equations (57) and (58) become, 
respectively. 




327T^CQr^ ^1 - cos 0^^ 



d^ dr 


3 


fcos^ B + cos^ 0^ 




‘(*|y>av = 


167T 




^1 - cos 0^ 



DISCUSSION OF EQUATIONS FOR FAR-FIELD IIMTENSITY; COMPARISON WITH 
LIGHTHILL'S THEORY AND RIBNER'S QUADRUPOLEMODa 

The first term on the right side in either equations (59) or (60) represents the self 
noise, and the other terms represent the shear noise (refs. 22 and 4). The factor 
^1 - cos 0^”^ in equation (60) represents the effect on the directivity pattern of the 
convection of the noise sources by the mean flow. The results of the present theory dif- 
fer from those of Lighthill's theory (as corrected by Ffowcs Williams in ref. 3) prin- 
cipally in that Lighthill obtains an exponent of -5 instead of -3 for this factor (refs. 2 
and 7). Ribner (ref. 5) refers to this term as the "convection factor. " Because of the 
large exponent, it is believed to be the most important factor in determining the direc- 

3 

tivity pattern of the jet noise. 

The term in brackets in equation (60) represent the actual generation term for the 
jet noise. It is composed of a self -noise part and a shear -noise part. The directivity 
pattern associated with this term is called by Ribner in reference 5 the "basic" direc- 


^This may not be true near the jet axis where refraction effects may be dominant. 
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tivity pattern of the jet noise. This pattern is a resvilt of the quadrupole nature of the 
noise source. 

In order to represent the "basic” directivity pattern given in equation (60), it was 
necessary to make some assumptions about the turbulence. These assumptions are the 
same as those made by Ribner (ref. 5) to deduce the basic directivity pattern implied by 
Lighthill's theory. Hence, it is appropriate to compare the basic directivity pattern ob- 
tained herein with that obtained by Ribner. The two results are almost identical. Thus, 
the self -noise term is the same in both cases and the directional dependence 

4 2 

cos 6 + cos 6 

2 

of the shear -noise term is also the same in both cases. The time derivative of the cor- 
relation volume multiplying (cos 0 + cos 0)/2 in Ribner 's model will reduce to the time 
derivative of the correlation volume 



obtained herein if the mean shear dU/dyg is assumed to be constant over this volume. 
Thus, the two shear -noise terms differ only by a small numerical factor. 

Ribner (ref. 5) attempted to carry his model further and estimated the ratio of the 
maximum shear noise to the self noise by making an assumption about the scalar F 
appearing in the turbulence correlation R^. Thus, (in the present notation) he assumed 

~9 ,r2t2/T2 

F = w^e"’^ ^ g(7) 

where L is a longitudinal correlation length and g(r) is a function of the time delay. 
With this assumption he shows that the ratio of the maximum shear noise to the self noise 
is proportional to 


9^g( t) 


1 

CD 

L 

7=0 

kg'(r)l 


1 

CO 

1 

7=0 


( 61 ) 


28 



And by assuming g(T) = e , he determines this ratio to be equal to 1/4. However, 
it is easy to show that upon choosing g(r) to be a decaying exponential times a low- 
degree polynomial, the coefficients of the polynomial can be adjusted to make the 
ratio (61) equal to zero or to any other value. Hence, it is not possible to use Ribner's 
argument with any confidence to estimate the ratio of the shear noise to the self noise. 
However, if it is assumed that the maximum shear noise is less than or equal to the self 
noise, the "basic" directivity associated with the noise sovtrces will be small compared 
with the directivity due to the convection factor ^1 - cos K will be seen subse- 

quently that this assumption leads to excellent agreement with the experimental data. 

COMPARISON WITH EXPERIMENT 

In this section we compare the predictions of the present theory with experimentally 
observed results on jet noise. The most striking comparison occurs when the theory is 
compared with the observed directivity data. 

The directivity patterns of the far -field intensity of noise from round jets have been 
measured by a number of investigators. We shall compare the present theoretical direc- 
tivity patterns with those found in the most recent experiments. In order to make this 
comparison, it is necessary to know the relation between the convection Mach number of 
the eddies M„ and the jet Mach number M,. Measurements of these Mach numbers are 
presented in references 15 and 26 to 28. There is general agreement among the various 
investigators. Their results show that the convection velocity is fairly constant across 
the mixing region and that in the center of the mixing region, where the most intense tur- 
bulence occxms, the relation = 0. 63 Mj holds. (As indicated previously, most of the 
turbulent energy lies in a narrow region centered along the center of the mixing region. ) 
We shall use this relation in comparing the theoretical results with experiment. 

As explained in the preceding section we shall neglect the effect of the shear noise on 
the directivity pattern, and therefore this pattern will be completely determined by the 
convection factor ^1 - M^ cos 0^ ^ . This quantity is compared in figures 1 to 4 with the 
recent experimental results of Olsen^ (ref. 29), Lush (ref. 11), and Krishnappa and 
Csanady (ref. 30) and with the older data of Howes (ref. 12). In all cases the level of the 
theoretical directivity curve is determined by pxitting it through the experimental data at 
90° to the jet axis, where the convection effect is zero. The curve corresponding to the 


“^These data are an improved version of the data in ref. 29 and will be presented as 
an AIAA paper. 
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present theory is shown as a solid line. For comparison purposes, the directivity factor 
predicted by Lighthill's theory is shown as a dashed line. In order to be as fair as pos- 
sible, the effects on the directivity of the shear noise are also neglected in the repre- 
sentation of Lighthill's theory. Including the shear noise would make the disagreement 
between the predictions of Lighthill's theory and the data greater in all cases shown. It 
can be seen from the figures that the present theory always agrees well with the data, 
and at least as well as Lighthill's theory at all jet Mach numbers. The agreement be- 
tween Lighthill's theory and the data is poor at the higher Mach numbers. Both theories 
and the data tend toward better agreement as the jet Mach number is reduced, primarily 
because the convection factor tends to unity as M^ approaches 0. 

CONCLUDING REMARKS 

A theory of aerodynamic noise has been developed which differs from Lighthill's 
theory primarily in the way the convection of the noise sources is treated. 

For high subsonic convection velocities, the present theory provides much better 
agreement with directivity measurements than does Lighthill's theory. For moderate 
and low subsonic velocities, both theories and the experiments tend toward agreement, 
as expected. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 3, 1972, 

501-04. 
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APPENDIX A 


EVALUATION OF 





In order to evaluate this integral, notice that 



for > 4 


for < 4 


where dr is the element of solid ai^le centered about ^ = 0 and ^ = |^|, I- 

Hence, if Sj^ is any sphere of radius R about ^ = 0, 
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APPENDIX B 


SYMBOLS 


A 

A 

B 

C 


c 


c 


o 


D 


D/Dt 

E 

e" 


F 

f 


g(r) 

1(^ 

I(x(^ 

A 

i 

L 

I 

Mj 

P 


P 

Po 

^nan 
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function in isotropic tensor (eq, (53)) 
vector potential 

function in isotropic tensor (eq. (53)) 
fimction in isotropic tensor 
specific heat at constant pressure' 
speed of sound 
average sound speed 

jet nozzle diameter, function in isotropic tensor 
0/at + V • V 

function in isotropic tensor 
viscous stress tensor 

arbitrary function of r, and t; fvinction of defined in eq. (54)) 
arbitrary function 

time dependence of correlation tensor 
far -field intensity 

intensity at x due to unit volume of turbulence at y 
spectral density of far -field intensity 

spectral density of far -field intensity of sound emitted from point y 

vinit vector in Xj -direction 

correlation length 

eddy correlation length 

convection Mach number 

jet -exit Mach number 

2 

Fourier transform of (p - Pq)/PqCq 
pressure 

average pressure in far field 

moving-axis, vector -potential correlation tensors 



R 

Ry(y,F, -t) 1 
Re 





U,U 


V V. 
» '^1 

V’ 


V, Vj 

W, Wj 


x,x. 


”^T T I T 

y,y ,y ,y ,z 


% 

6/6y. 

^ijk 

’Ji 


gas constant; radius of sphere about ^ = 0 

moving -axis velocity correlation tensors 

moving-axis, turbulence velocity correlation tensor 
real part 

fixed -axis turbulent correlation tensor 
distance from source point to field point, |x - ^ 
component of r ; Xj^ - y^ 
entropy 

large time duration 
Fourier transform of t- 
time 

mean velocity 

convection velocity 

total fluid velocity 

fluctuating part of velocity 

acoustic part of fluctuating velocity 

incompressible velocity 

turbulent part of fluctuating velocity 

Cartesian coordinates; location of observation point 

location of source point 

(0/at + u 0/3xj)n 

ratio of specific heats 
Krone cker delta 
partial derivative with r fixed 
permutation tensor 
z - y 

-nt 

z - y 
-y” 

component of rj 



0 temperatvire 

0 angle between mean flow direction and observation direction 

1 

7^^-U^(y)T 

^2) ;j(2) _ 

component of ^ 
n y“^ In (p/Pq) 

p density 

p' fluctuating part of density 

Pq average density 

T time delay 

Tj characteristic decay time of moving -axis correlation 
Tjj turbulent stress tensor 

q> azimuthal angle 

oj(l - M cos e) 

CO circular frequency 

coj characteristic frequency of second-order correlation 

Superscripts: 

' quantity evaluated at ^ and t + t 

” quantity evaluated at y” and t 

quantity evaluated at and t + t 
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Intensity, dB 






Figure 3. - Experimental directivity 
data of reference 30. Jet nozzle 
diameter, 2. 54 centimeters 
(1 in. ). 



(a) Jet nozzle diameter, 12. 7 centimeters (5 in. ). (b) Jet nozzle diameter, 7. 62 centimeters (3 in. ). 

Figure 4. - Experimental directivity data of reference 12. 
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